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I. INTRODUCTION 

The proton isovector-axial coupling and quark mo¬ 
mentum fraction {x)u-d are important benchmarks to 
check whether the systematic uncertainties of lattice 
QCD simulation, such as finite lattice spacing, hnite vol¬ 
ume, and chiral extrapolation, are under control, by a 
correct reproduction of the corresponding experimental 
results. Since the noisy disconnected insertion contribu¬ 
tion to the isovector part of the nuclear matrix element is 
canceled between two degenerate flavors, the values are 
obtained solely from the connected insertion and thus are 
relatively cheaper to compute with high precision to be 
considered as benchmarks. 

Most attempts have resulted in values ~10% below the 
experimental number for the axial-vector coupling m, 
while a few claim that their results could be consistent 
with experiment [mg. For the quark momentum frac¬ 
tion {x)u-d, overestimation by ~20 - 30% is common in 
most of the calculations o [3 [HHISl except [5]. 

Recently, attention has been paid to lattice QCD cal¬ 
culation of the isovector scalar matrix element gg in the 
proton la mi nil ng due to its role in constraining pos¬ 
sible scalar interactions at the TeV scale [15] . 

In this work, we calculate the isovector matrix elements 
of the nucleon for the axial-vector and scalar couplings 
and the quark momentum fraction with the valence over¬ 
lap fermion on 2 -|- 1 flavor domain-wall fermion (DWF) 
configurations [20) . Compared to simulations with other 
actions, the overlap fermion provides the best control of 


the systematic errors since it is free of explicit chiral sym¬ 
metry breaking and gives small 0 {a^) errors, whereas the 
numerical work is more costly. 

In order to improve SNR, the 8-grid smeared Z 3 noise 
source with low-mode substitution (LMS) [55H57] has 
been applied to the hadron two point correlator on the 
24^ X 64 lattice [3g which improves the error of the nu¬ 
cleon mass of a point source by a factor of 7 and that 
of the 8-grid source without smearing by a factor of 2.5. 
In this work, we use a stochastic sandwich contraction 
method to remove the need of multiple inversions in the 
sink-sequential approach and use the current-sequential 
method for the low modes in the propagator between the 
current and the sink. This is an extension of the noise 
grid smeared source with LMS to the three point func¬ 
tion. Such a many-to-all correlator with LMS is useful 
when the low-eigenmode contributions are important in 
the relevant time windows where the physical quantities 
are extracted. 


The structure of the rest of the paper is organized as 
follows. The LMS technique with noise grid source for 
the non-zero momentum case of the two point correla¬ 
tion function is provided in Sec. [Tlj Sec. O discusses the 
possibility of applying LMS on all the four quark propa¬ 
gators in the proton three-point function. The numerical 
details are provided in Sec. [TV] In Sec. jV] the results of 
isovector matrix elements of the nucleon for the axial- 
vector g\, the scalar coupling gg and the quark momen¬ 
tum fraction {x)u-d are provided. A short summary and 


outlook are presented in Sec. VI 
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II. LOW MODE SUBSTITUTION WITH MIXED 
MOMENTUM GRID SOURCE 


with Cc as the upper bound of the modulus of the eigen¬ 
values. 


Let’s first consider the nucleon two-point function 
(2pt) with the interpolation field of the nucleon |29) . 


Xo.{x) = 

Xa'{x) = -e (x), (1) 


where C = = 727475 in the Pauli-Sakurai gamma- 

matrix convention, used throughout this work. There 
are two kinds of the Wick contractions so the 2pt of the 
nucleon can be constructed in terms of the point-to-point 
quark propagator S as 

Tr (y, x)) Tr (y, (y, x))) 

- (e“''^e“''>''='Tr {y, {y, x)S'(“)“' {y, x) j) 

_ ^^abc^ab'c 

Tr (y, x)) Tr (y, (y, x)) 

+Tr (y, (y, x)5(“)“' (y, x))) (2) 


where «S is defined as ((75(7“^)^ and T is the projection 
operator for the nucleon polarization. 

The quark propagator S in the above equation is the 
inverse of the operator {Dc + m) [301 El] , where is 
defined in terms of the overlap operator and is chiral, 
i.e. {Dc, 75 } = 0 [52|. The details will be discussed in 
Sec. 1^ As in Ref. [3T1 [35] , we use the low lying eigen¬ 
values and eigenvectors of the overlap fermion, Ai and 
|z), satisfying Dc\i) = Ai|j) to speed up the inversion and 
separate the propagator into its low-mode and high-mode 
parts. 


5 'l( 2 /,x) 


SH{y,x) 


E 

|Ai|<e, 


1 




S{y,x)- ^ 


1 


|A.|<e 


\i + m 


l*)y(*U) (3) 


The idea of using the Z 3 noise grid source is to tie 
the sources of the three quark propagators stochastically 
to each point (or a smeared point) on the grid so that 
one can have a multi-to-all correlator from one inversion. 
LMS for the quark propagator with Z 3 noise grid source 
(PropNG), be it point-grid (PG) [3T] or smeared grid 
(SG) [35], has been used to improve the SNR for the nu¬ 
cleon correlator with significant success. This technique 
removes the gauge non-invariant contributions of the low¬ 
mode contributions (defined below) from the cases in 
which three propagators are from different source sites, 
and restores the benefit of using PropNG. 

To construct the nucleon correlation function with 
LMS, PropNG SNciv) should be split into its high-mode 
and low-mode pieces 


SNciy) = '^^ix)S{y,x) 

xGG 

= ^^^(y)-f ^ 6»(x)S''^(y,x), (4) 

xGG 


with S§Q{y) = random Z 3 

phases 0 {x) S {1, e's’’’, for each point on a grid 

G. 


As in Ref. [28], we can expand the nucleon correlation 
function C'(y,x; P; S'^g, S^q, 5^^) with the decomposi¬ 
tion in Eq. @ (ignoring the indices for the sink position 
y and the projection matrix P), 


/^LMS f o Q c \ 

^ [^NG^ ^NG, ^NG) = 

= C{S^c, SSg. S§g) + E C{0ix)S^ix), e{x)S\x), e{x)S\x)) 

x^G 

+C{ Y. 0(x)5^(x), S^a) + C{S^g, Y S^g) + C{S^g^ S§g, Y 

x£G xGG xGG 

+ Y C{9{x)S^{x),0{x)S^{x),S^g) + Y C{0{x)S^{x),S^G,&{^)S^{^)) + E C{S^G,&i^)S^i^),0{x)S^ix)) 

x£G xGG xGG 

= CkeriS^G, Y Hx)SHx)) + Y Cker{0{x)S^ (x), S^g) 

xGG xGG 


( 5 ) 
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where 

Cker{Sl,S2) = C{Si,Si,Si)+C{S2,Si,Si) 

+C(5i,52,5i) + C(5i,^i,52). (6) 

The nucleon correlator with LMS here can be obtained 
from the one in Ref. [28j with just one more step. The 
low-mode propagator decomposed 

into several terms as in the very last term in the RHS of 
Eq. to improve the SNR. 

After the noise averaging, the nucleon correlation func¬ 
tion with PropNG should be a stochastic estimate of the 
sum of nucleon correlators from each of the grid points, 
i.e. 


'^Cgridiy) = Y^'^C{y,Wi), (7) 

V i V 

where the grid points Wi are 

Wi e {xo + rrixA^, yo + iriyAy, zq + m^A^). (8) 

with rnx,y,z = (0,1, • • • , Ls/A^^y^z) modulo the periodic 
boundary condition in the spatial directions. In this grid 
pattern, in addition to the zero momentum mode (0,0,0), 
one can obtain non-zero momentum modes from the nu¬ 
cleon correlation function with PropNG. For example, for 
the PropNG with a regular (A^, = Ay = = Ls/m) 

grid, the momentum mode p = (inim, ±n 2 m, inam) 
(^^ 1 , 2,3 are integers) can be obtained. In this case, there 
is a phase factor which needs to be taken into account 
when the origin wq = (xg, yoi^o) is changed from config¬ 
uration to configuration, 


E 








= e 


o.pEEC'(y,u;.)< 




(9) 


I y 


The exponential term in the second line with the expo¬ 
nent proportional to Wi — wg does not contribute, since 
all components of the latter are proportional to Lg/m 
and, as a result, the exponent is a multiple of 27r. 

In order to obtain the other momentum modes, propa¬ 
gators with noise grid non-zero momentum source (Prop- 
NGM) are required. To cover a range of modes and 
minimize the effect of the rotation symmetry breaking 
due to the finite lattice spacing and volume, three kinds 
of PropNGM 

SpAy) = 

i 

SpAy) = 

i 

^Psiy) = ^0{wi)S{y,Wi)e^^''^°’°’^^ ( 10 ) 


and related inversions are required for the proton case. 
It is trivial to confirm that one can obtain a momentum 
mode like (1,1,0) from the contraction C{Sp-^, Sp^, Snc), 
and (1,1)1) from C{Sp^, Sp^, SpA- 

To reduce the cost, we can combine these three kinds 
of PropNGM together as the mixed PropNGM, 

Sp = Sp^ Sp2 Sp^ 

= 9{wi)S{y, 

i 

with the origin of the grid wg = {xg, yg, zg) to be selected 
randomly for each configuration. 
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FIG. 1: The plot shows the relative error of 2pt as a function 
of the momentum squared at t=8 in lattice units. The data 
points of the smeared grid cases have been shifted a bit on 
the abscissa to make it easier to distinguish them. The SNR 
of the case with the noise smeared grid source (red squares) 
and LMS applied is better than the one with smeared point 
source (blue dots), while the one with the noise smeared grid 
source but no LMS (black triangles) is even worse than the 
one with smeared point source. 


Fig.g shows the SNR of the proton effective mass at 
the unitary point where the pion mass due to the valence 
quark is the same as that from the sea, on the ensemble 
of which details will be addressed in Sec. [IVl When LMS 
is applied, the SNR of the 2pt with the noise smeared 
grid source propagators (PropNG and mixed PropNGM, 
Ax = Ay = Az = Lsl2) is 2.3 times smaller than that 
of the of the smeared point source at = 0. This is a 
gain of 5.3 in statistics which is very good considering 
that the maximum possible gain is 8 for the ideal case 
where the independent nucleon propagators emerge from 
each of the 8 smeared grid points. On the other hand, 
if we don’t use LMS, the SNR of 2pt with grid source 
is worse than the smeared point source, even though the 
latter has only 1/8 of the statistics of the former. This 
is understood as due to the fact that the Parisi-Lepage 
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estimate of the SNR for the nucleon is modified to 

CN{t,P= 0) _ -(mN-3/2m„)t 

^Nit) ^ ’ 

where N is the product of the number of noise and 
the number of gauge configurations and V3 is the three- 
volume of the noise with its support on a time slice. In 
our case, V3 = 8. It is this extra factor of which 
makes the SRN of the 2pt from the noise smeared grid 
source without LMS worse than that of the smeared point 
source. When LMS is employed, the situation is reversed 
and one gains a statistical factor almost as large as the 


number of the grid points. Thus, it is essential to have 
LMS when the noise grid source is used for the nucleon. 


III. LMS OF THE CONNECTED THREE-POINT 
CORRELATOR 

Generally, a nucleon three point function (3pt), from 
X to y, with a current (with current 

operator O such as 7^, ^iDj, etc.) inserted at z, includes 
four kinds of Wick contractions. 


(7“ (y, x; L; ) 


l^^abc^a-b'c'^^ (^s{n)ad^y^ (z, x)) Tr (y, x)S'(“)“' (y, x))) 

^^^abc^a'b'c'rj.^^ (z, x)S^‘^'^^^' (y, (y, x))) 

^(^^abc^a'b'c'j,^ (rs'(“)““'(y,x)) Tr (^(‘')“'(y,x)S'(“)‘=‘^(y,z)C>(z)S'(“)‘^^'(z,x))) 

-f (e“'’^e“''’'^'Tr (y, x)S^‘‘'>^^' (y, x)S'(“)‘='^(y, z)C>(z)S'(“)^'=' (z, x))) (13) 


and can be expressed in terms of the 2pt correlation func¬ 
tion C'(y, x; T; S'(“^) defined in Eq. ([^, 

C'^(y,x;r;5(“),S'(“),S'(‘^\S'(“)) 

= C'(y,x;r;5(“\S'(''\S'(“)) 

+ C'(y,x;r;^(“),^W,5(“)), (14) 


where S{0,zo\y,x) = J2gS{y,z)0{z)S{z,x) is the cur¬ 
rent inserted propagator (PropCI). Similarly, the 3pt 
with a current of d quark can be expressed as 


PropCI (from 0 to t 2 through 
z(z,ti) 



C^(y,x;r;5(‘^),S'(“),S'W,S'(“)) 

= C'(y,x;r;S'(“),5(‘'),S'(“)). (15) 


Fig. [2] shows PropCI as the product of the propagators 
in the shadowed region. 

Supposing = S', Eq. 

into the contraction of PropCI S and the remaining parts 
denoted as A„_d(r, Si, S2), 
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can be rewritten 


FIG. 2: The quark diagram of the proton correlation function 
with the connected insertion, from x to y, with an insertion 
at z. The product of the propagators in the shadowed region 
is the current inserted propagator, S. The propagator from 
the current 2 to the sink y is decomposed into its low- and 
high-mode contributions (Sl and Sh respectively) for further 
SNR/cost improvement from the advanced technique in the 
latter discussion. See Sec. IIIIBI for more details. 


C 3 “(r;S,S,S,S) = (Tr(SX,(r,S,S))), 
C 3 "(r;S,S,S,S) = (Tr(SA,(r,S,S))), (16) 
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with 


X aa' C c \ ^abc^a'b'c' 

u *31,02; = e e 

(rTr[^“'5“'] +^“'s'r'r 
+Tr[rs'“']^f + 

X bb' /T' C C ^ ^abc^a'b'c 

d , *31, *32j = e e 

(Tr[rS'r']C-i(S'2“')^C' 

+C-\Si°^°^'TSffC) (17) 

Based on the above definition, a typical 3pt correlation 
function for a point source on the t = 0 time slice, when 
summed over the spatial indices of y and z becomes 

C 3 (t 2 ,ti) = ^ (Tr[5(C>,ti;y,t2,0,0) 
y 

^«.d(y,i2,0,0;r,^,^)]). (18) 


A. Sink-sequential method and Stochastic 
sandwich method 

The typical problem of the connected 3pt is calculating 
the propagator from the current to the sink S{y, t 2 ,z,ti). 
On the surface, it is an all-to-all propagator which would 
be beyond the ability of the standard lattice inversion 
operation. 

However, when the sink time O is fixed, the se¬ 
quential source method [3S1 El] could be used, with 
75 ^(y, 0 , 0 , 0)75 as the source of the matrix inversion, 
to construct 

Sseq{Xu,d;z,ti,t2,0,0) = y^5'(z, 0,^,0) 

y 

^ 2 , 0 , 0 ) 75 . (19) 


many Np are required for nucleon form factors with mo¬ 
mentum transfer (hundreds are needed for |p| < 3 with 
high statistics), the cost can be staggering. 

A stochastic method [57H55] (referred to as the 
stochastic sandwich method (SSM) in this work) is intro¬ 
duced to reduce the cost of the sequential method when 
many sequential inversions are required. It entails in¬ 
serting a noise estimate of the delta function i5(yi, 1 / 2 ) at 
t = t2, 


, iVnoi 

X! Tr[6»l’^^S'(yi,t2,2,ti)0(^,H)*S'(z,ti,0,0) 

^(0,0,y2,^2)4^^] J^-^C 3 {t 2 ,h, 0 ), (21) 

where Nnoi is the number of the noises and the noise 9 
satisfies 






y 

yi 2/2 


iVnoi-)-00 


( 22 ) 


i=l 


In other words, it uses the noise estimate of the all-to-all 
propagator, 

Siyi,t 2 ,z,ti) = y (23) 


with 

Snhiz,ti,t2) = y S'(z,^l,yl,^2)6'^*^^ (24) 

VI 

instead of the original S(y,t 2 , z,ti), to avoid the expen¬ 
sive calculation to construct the sink-sequential propaga¬ 
tor with inversion of 2 x 4 x Np sources. 


B. Stochastic sandwich method (SSM) with LMS 


Then, one can contract Sseq with the standard quark 
propagator from t = 0 to H to construct the 3pt correla¬ 
tor, 

C 3 {t 2 ,ti, 0 ) = yTr[ 75 y^(A„,d, 2, H,i2,0,0)75 

z,i 

Oiz,h)Siz,h,6,0)], (20) 

taking the advantage of the relation 75 <S'(z, y)i 75 = 
S{y,z). 

The disadvantage of the sequential method is that it 
has to calculate the sink-sequential propagator repeat¬ 
edly when X is changed for any reason, such as for: 
different momentum, different quark flavor or mass, or 
different polarization projection of the baryon. This is 
expensive when many momenta are needed. 

The number of inversions required in the sink- 
sequential method is 2 x 4 x Np where the 2 is for the 
u and d flavors in the nucleon, 4 is for the polarization, 
and Np is the number of momentum projections. When 


SSM avoids the cost of the repeated inversion for many 
different sequential sources, but it still requires multiple 
inversions for several noises, before the SNR can reach its 
upper limit - that of the sequential method. In this work, 
the basic idea is to improve the SNR of the 3pt correlator 
of SSM using the low lying eigenvectors of Dc to construct 
the long distance part of the all-to-all S{y,t 2 , z,ti) {Sl 
in Fig. the single line from the current to the sink), 
and using the noise many-to-all propagator to estimate 
the remaining high frequency part of S{y,t 2 , z,ti) {Sh 
in Fig. the double line from the current to the sink). 
Thus, the propagator with LMS is written as 

= y 6l^^^75(*S'io)^(5',ll,i2))^75 

i 

+ y ' ^ Vt[yM)vl{z,ti). (25) 

~r 'fTl 

i 

where Xi and Vi are the low-lying eigenvalues and the 
corresponding eigenvectors of Dc- In other words, it is 
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a technique to apply LMS to the sequential propagator 
Sseq{Xu,d': ^ 1 , ^ 2 ,0,0) (LMSg). It is expected to reduce 
the number of the noise propagators needed to reach the 
upper limit of SNR. 

When LMSs in Eq. (25) is applied to the PropCI in 
Eq. (14), S comes from t = 0 to t = t 2 through t = ti 


as shown in the shadowed area in Fig. 


S^^^^iO,h;y,t2,h,6,0) = 

= {yi,t 2 , z,ti)O{z,ti)S{z,ti,Q,0) 

Z 

z i 

+0''"Ky, t2) ^ lb{S^nlf{z, ti, t2))^75) 

z,i 

O{z,ti)Siz,h,6,0), (26) 


Then one can construct 3pt with LMS by constructing 
the standard 2pt repeatedly (the projection matrix E is 
suppressed for clarity), 


C. 


LMS,u 


^LMS.d(^^^) 


iS,S) = ClAS^,S\Sf}c,SkG,SSG.SkG) + E ClAS^^^\S^,0i^)S\S§G,0i^)S\S§G) + 

xGG 

ClAS^^G,SkG.S^G,SkG,S^,S^) + E ClMx)S^SSG,di^)S^,SSG,S^^^\S^) 

xGG 

CLriS^G^S^GJ'',S^,S^G^SkG) + E CtM^)S\S^G^S^<^^\S^,0{x)S^,S^G) (27) 


xGG 


r 


where 

S^G = E y{^)S^{x), and 

xGG 

Ct,iXi,X2,Yi,Y2,Zi,Z2) = CiXi,Yi,Zi) 

+C{X2^ Yi, Zi) + C{Xi,Y 2, Zi) + C{Xi,Yi^ Z 2 ) 

(28) 

and and are the high- and low-mode parts of 

gLMSs in Eq. (1^. 

This is the stochastic sandwich method with LMS 
which uses the low eigenmodes for the propagator from 
the current to the sink in PropCI, with current in¬ 

sertion and the high modes for the same which originates 
from the sink time slice. The construction of the PropCI 
with low modes needs to be done for each current and 
momentum transfer and ^2 (if desired). In contrast, the 
current-sequential method will need to do an inversion 
for each current, momentum transfer, and ti separately. 

To account for the amount of numerical work for differ¬ 
ent approaches to the 3pt Cl correlators, we note the the 
traditional sink-sequential method entails 2x4xiVp inver¬ 
sions at a fixed sink time slice ^ 2 , where the 2 and 4 refer 
to the separate sources X in Eq. (17) labeled with u and 


there are Nnoi inversions of the Nnoi noise vectors at the 
sink time t 2 . How many Nnoi is needed for acceptable 
SNR depends on the observable. For the SSM with LMS, 
besides the noise propagator Snoi with inversion, 

there is an over head for the low-mode portion of PropCI 
(^gLMSs jjj Eq. (26l). It includes N times the low-mode 
contributions from N smeared grid source plus one high¬ 
mode contribution for the propagator from the source to 
the current (S'^g). Each needs to be folded with the 
current for different momentum transfer q. Therefore 
the overhead is e x (TV -|- 1) x Ncu x Nq where Ncu/Nq 
is the number of currents/momentum transfer, and e is 
the fraction of inversion time for constructing the low¬ 
mode portion of for each current and momentum 

transfer. We list the cost for the sink and current parts 
of the 3pt function in units of quark inversion in Table |T] 
for future reference. To evaluate the efficacy among the 
three methods, one needs to compare costs in the table to 
reach the same precision for a given observable. For the 
case of SSM with LMS, there is an additional gain from 
the noise grid source with LMS as discussed in Sec. [IT] 
which needs to be taken into account. 


d flavors and polarization directions (unpolarized and po¬ 
larization in 3 spatial directions). Np is the number of 
sink momenta for the nucleon. For SSM without LMS, 
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TABLE I: The cost for the sink and current parts of the 
3pt function in units of quark inversion is listed for the sink- 
sequential method (Sequential), stochastic sandwich method 
(SSM), and SSM with LMS. Np is the number of sink nu¬ 
cleon momenta, Nnoi is the number of noise in SSM. is 
the number of noise in SSM with LMS, and Ncu/Nq is the 
number of currents/momentum transfer in the construction 
of of the low-mode part of PropCI. e is the fraction of inver¬ 
sion time for constructing the low-mode portion of PropCI 
for each current and momentum transfer and Np momenta 
(~ 0.02 on the ensemble used in this work). 


Sequential 

SSM 

SSM-kLMSs 

00 

^noi 

N^ciMN + l)NcuNq 


IV. NUMERICAL DETAILS 

In this work, we use the valence overlap fermion on 
2-1-1 flavor domain-wall fermion (DWF) configurations 
m to carry out the calculation m- 

The lattice we use has a size 24^ x 64 with lattice 
spacing a~^ = 1.75(4) GeV set by vq at the chiral 
and continuum limits [42) . The light sea u/d quark 
mass mia = 0.005 corresponds to ~ 330 MeV. We 
have calculated the isovector matrix elements of the nu¬ 
cleon for the axial-vector and scalar couplings and the 
quark momentum fraction at 6 valence quark mass pa¬ 
rameters which correspond to the renormalized masses 
= TO^®(2GeV) ranging from 13 to 32 MeV after the 
non-perturbative renormalization procedure in Ref. |41j . 
They correspond to the pion mass in the range of 250- 
400 MeV. In order to enhance the signal-to-noise ratio 
in the calculation of three-point functions, we use two 
smeared noise 12-12-12 grid sources at ti = 0 and 32 
(one is PropNG and the one is PropNGM) [55] and two 
noise 2 - 2-2 grid point sources at positions tf which are 
8 , 10, and 12 time-slices away from the sources on 203 
configurations. 

The effective overlap operator Dc is chiral, i.e. 
{Dc, 75 } = 0 [55], and is expressed in terms of the overlap 
operator Dov as 

Dc = ^ with Dov = l + 75 e( 75 -Dw(p)), (29) 

where e is the matrix sign function and is the Wilson 
Dirac operator with a negative mass characterized by the 
parameter p = 4 — 1/2k for Kc < k < 0.25. We set k= 0.2 
which corresponds to p = 1.5. 

Compared to the earlier implementation of the over¬ 
lap operator | 21 j , the current implementation further im¬ 
proves the performance of data exchange on different 
nodes of the cluster and uses the polynomial approxi¬ 
mation for the overlap operator instead of the rational 
approximation, and has achieved better scaling and fur¬ 
ther speed up of the calculation by a factor of two on 
average [55] . 


The number of D^’s low mode eigenvectors used for 
the deflation of the overlap operator inversion and LMS, 
on this 24^ x 64 lattice, is 200 pairs plus the zero modes, 
and the upper bound of the absolute value of the eigen¬ 
values is 0.154 which is over two times larger than the 
dimensionless strange quark mass. 

We check the efficacy of the sequential low-mode sub¬ 
stitution (LMSs) in the PropCI by examining the 3pt 
functions for the isovector axial and scalar currents. We 
plot the ratio of 3pt-to-2pt correlators as a function of 
the current insertion time ti in Fig. [^ where the sink 
time t 2 is 10. The blue dots and black triangles show the 
contributions where the current-to-sink part of PropCI is 
from the low modes and the noise-estimated high modes 
respectively. Notice that the contribution from the low 
modes is much larger than that of the high modes when 
the current time slice is farther away from the sink (i.e. 
closer to the source with small ti) for both the axial and 
scalar cases, which reflects the fact that the low modes 
dominate the long-distance behavior of the PropCI be¬ 
tween ti to t 2 . When the current is closer to the sink with 
larger ti , we see that the high modes dominate for the ax¬ 
ial case which shows that the high modes are important 
and dominate the short distance behavior of the propa¬ 
gator. However, the high-mode contribution is still small 
for the scalar current case when ti is close to the sink 
which shows that the high-mode contribution is small 
for the 3pt function for the scalar current. 

The red squares are the sum of the low- and high¬ 
mode contributions from the present hybrid scheme. We 
have also calculated the 3pt function without LMSs for 
the PropCI, but instead use only the noise propagator 
as the full propagator from ti to ^ 2 - These are shown as 
the green triangles in Fig. [^ These correspond to the 
stochastic method introduced by the QCDSF Collabora¬ 
tion [37| [38] and the Cyprus group [39]. Since our LMSs 
replaces the long distance part of the current-to-sink part 
of PropCI with an exact all-to-all one, the larger its con¬ 
tribution the larger the improvement. As in Fig. [^ the 
blue dotd contribute over 80% in the gg case and so the 
improvement of LMSs is larger than in the g\ case. The 
error bars of SSM at the time slices C = 2 — 6 turn out 
to be a factor ^ 2 for g\ (^ 4 for gg) larger than that of 
using LMSs in the present approach. 

The fact that the error of g\/g% in our approach is 
smaller than that of SSM with 2 noises by a factor of 
~ 2/4 shows that it would take 8/32 noise inversions for 
SSM to have the same error as the present method with 
LMS. To compare the cost of SSM + LMS, we should 
take its overhead into account. On the present lattice, the 
percentage of inversion time for low-mode construction is 
e = 0.02. Therefore, the overhead e(iV -|- l)NcuNq = 0.72 
for = 8 (smeared grid source), Ncu = 4 to account for 
the scalar current and Ai for 3 spatial directions and 
Nq = 1. Together with = 2, the cost is 2.72 inver¬ 
sions. This means that, to reach the same error, it would 
take SSM 2.9 and 11.8 times more inversions than SSM 
with LMS for g\ and g^g respectively. Furthermore, the 
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FIG. 3: The 3pt-to-2pt ratio with LMSs (red squares) vs. the 
one without it (green inverted triangles). The source/sink is 
located at 0/10, and the current dependence for the matrix 
element with the current-to-sink part of PropCV including 
just the low- or high-mode parts are plotted as the blue dots 
or black triangles respectively. The upper panel is for the 
axial-vector current case and the lower panel is for the scalar 
case. Notice that the contribution of the low-mode part is 
larger when the current time slice is farther away from the 
sink. 


smeared grid source with LMS has improved the statis¬ 
tics by a factor of 5.3 for = 8 for the 2pt function. 
This additional factor of improvement is also expected 
for the 3pt function. 

To compare with the sink-sequential method, we as¬ 
sume that our results have reached the SNR of that of 
the sink-sequential method. This is consistent with the 
fact that in the range ti = 2 — 6 where the observables are 
fitted, the PropCI are dominated by the low-mode con¬ 
tributions, particularly for g^. In this case, the cost of 
sink-sequential takes 16 inversions. Here, we have taken 
Np = 2 to include the {x)u-d calculation in addition to 
g\ and g^g. For the overhead in SSM -|- LMS, the number 
of currents needed is N^u = 6 for these three quantities 


FIG. 4: The vector renormalization constant in the 
rest/moving frame at the unitary point, as a function of the 
momentum squared in lattice unit. The p^ =0,4 involve 
PropNG only and the other cases involve PropNGM also. The 
former case gains more from LMS (black squares vs. red dots). 
The results obtained using the 8-point grid source without 
LMS applied are very noisy (blue triangles). 


and the overhead is e(7V -|- l)NcuNq = 1.08. Therefore, 
besides the improvement from use of the grid source, the 
present method would be 16/3.08 = 5.2 times more effi¬ 
cient than the sink-sequential method for the calculation 
of the three quantities. Note that the cost of the sink- 
sequential method has additional factors that need to be 
taken into account, such as iV^ass for different masses, 
and also N when the necessary LMS is applied on the 
source of the sink-sequential propagator (as in Eq. 171, 
so SSM is much cheaper than the sink-sequential method. 

When the physical volume is increased, while keeping 
the lattice spacing unchanged, and with a noise vector 
covering the entire spatial volume of the sink time slice, 
we expect that the region essentially contributing to 3pt 
will not change, while the remaining region contributes 
only to the noise. Such a simple argument hints that 
the noise required to reach the same SNR is proportional 
to volume and we have confirmed it explicitly on the 
48^ X 96 lattice with similar lattice spacing [?D]. At the 
same time, the number of low modes will be proportional 
to volume if we want to reach the same upper bound of 
the eigenvalues, so the SSM with LMS will not lose its 
efficiency as compared to SSM without LMS, when the 
volume is larger. But, since the number of inversions is 
fixed in the standard sequential method, the SSM with 
and without LMS will lose their comparative efficiencies 
when the volume is very large. 

Another issue we need to check is the effect of LMS in 
the 3pt case. For the 3pt function, we check, for exam¬ 
ple, the vector charge renormalization constant from the 
forward matrix element at the unitary point for several 
nucleon momenta. For = 0 and 4, only the propagator 
PropNG is involved, while the other cases involve Prop- 
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NGM also. In the former cases, we find that the smeared 
grid source with LMS improves the SNR by a factor of 
2.0 compared to that with a smeared point source with¬ 
out LMS, slightly smaller than what we found with the 
2pt function as discussed in Sec. whereas, the gain is 
only 1.4 for the other where the PropNGM is involved. 
We shall look into the possibility of improving the SNR 
further when PropNGM is involved. 


V. RESULTS 


three values of t 2 so that we can fit all four terms in 
Eq. 

In order to check the ^2 dependence of the plateau, 
three sets of propagators with two noise-grid point 
sources each at positions ^2 = 8,10 and 12 time-slices 
away from the nucleon source are generated, and all the 
ti dependence of these three cases are plotted together 
for comparison in Fig. for the vector current case. The 
sink-source separation dependence seems to be mild here, 
but in general the minimum separation required by other 
quantities can be different. 


A standard 3pt/2pt ratio in the forward matrix ele¬ 
ment case is 

R(t2,ti,0) = C3{t2,ti,0)/C{t2,0) 


(fc) y{k) 




(xri'^lxD 


( 0 )\ 


.(1) 






(1) 


(0) 




-b 


v(l) v(l) 

-f5)^((x['Vlx['^) - (x[°Vlx[°^))e-^^*^ 


+..., (30) 

where and are the energy and the overlap of the 
interpolation field of the ith state and AE = — E^°\ 

For t2 ^ ^ 0, the contributions from all the terms 


in the right hand of Eq. (30) except the first term van¬ 
ish, and then one can use Eq. (30) to obtain the matrix 
element. 

When t 2 is fixed, one may fit the first term and the 
combined second and third terms around ti = t 2 /2 to in¬ 
clude the effect of the ground state to first excited state 
transition in the right hand side of Eq. (30) which is ti 


dependent. But since the fourth term in the right hand 
side of Eq. (30), which is the difference of the matrix 


element in the ground state and the first excited state, 
is independent of ti just like the first term, one will not 
be able disentangle them and, as a result, a systematic 
error may be induced by its contribution which is sup¬ 
pressed by feeling for the size 

of the correction, let us suppose that the first excited 
state matrix element (Xf^^I'^lxP^) is 30% different from 

the ground state matrix element (Xf°^I'^IXi^^): and the 
mass difference of the first excited state and the ground 
state is about 500 MeV. Then the correction from such 
a effect with ^2=8, 10 and 12 (with the nucleon source 
set at to = 0) is about 3%, 2% and 1% respectively. To 
assess this error, we shall calculate the 3pt function at 



0.8 - ^^^^^^ 

-6 -4 -2 0 2 4 6 

H2/2 

FIG. 5: The nucleon sink-source separation dependence of 
the matrix element with the vector charge for u -|- d in the 
connected insertion. Obviously, the larger t 2 , the worse the 
signal. The data points marked with the black squares (t 2 = 8 ), 
the blue dots (t 2 = 10 ) and the red triangles (t 2 = 12 ) are con¬ 
sistent. 


To check the separation effect quantitatively, we ap¬ 
plied three kinds of fits to deduce the results: 

The first method is to fit the ratio as a function of ti 
and t 2 , 

RE\t2,ti) = Co 

-bC2e-^'"‘i -b C3e-^™‘== (31) 

with Co,1,2,3 and Am as free parameters. Co is the 
ground state matrix element we want. Since the C de¬ 
pendence of R{t 2 ,ti) is mild in some of the quantities like 
gv and qa, we take Am as a common parameter for all 
the quantities. This is what we mark as “2-state” in the 
following discussion. 

In this work, we use the smeared source and the point 
sink, so the excited-state contaminations are different in 
the smaller and larger t ends. If the smeared source 
makes the contaminations in the smaller t end small, or 
has a different sign compared to that in the larger t end, 
the position of the plateau will be harder to determine, 
as in the case of (Fig. and g\ (Fig. [^. Ap¬ 
plying the “2-state” fit on such a quantity is not stable 
and provides large uncertainties (and/or large /d-o.f.) 
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on the results. In this work, we constrain the mass dif¬ 
ference Am to be the same for the different matrix ele¬ 
ments with the same quark mass value, and apply a cor¬ 
related joint 2-state fit. To suppress the contamination 
from the excited state, we excluded the data points with 
ti — 0,t2 — 1 and <2- One more data point at the larger t 
end is excluded since the excited-state contamination is 
larger there. Despite this, the fit is still not very good. 
Taking the unitary point as an example, the x^/d.o./. 
with ~70 degrees of freedom is 1.45, the corresponding 
p-value is just 0.008. In addition, this method requires 
a joint fit with several quantities and is not suitable for 
the analysis of a single quantity without the information 
of the other quantities. 

The second method is the sum method [321111] which is 
used in the disconnected insertion case, wherein a sum is 
taken over all the 3pt/2pt ratios in Eq. (30) with different 
ti, 


SR{t 2 ,ti, 0 ) — i?(t2,ti,0) 

0<tl<t2 

= {h - 

dl) 


^ — Am . 7 ( 1 ) , . 


7(1) 7(1) 

+(^2 - 1 ) J(o) - (x[°Vlx®)) 

Zf Z. 

-Amt 2 1 


(32) 


When t 2 is large, we can use the linear function of t 2 
(ignoring the correction) 

SR^^\t2,h) = t2Co + C[ (33) 

to fit our summed ratio with 3 different separations, and 
obtain the slope as the ground state matrix element. This 
method will be marked as “sum” in the following discus¬ 
sion. 

We found that the “sum” fit can obtain a ^ jd.o.f. 
smaller than one, for all the quantities. But this fit just 
has one degree of freedom. Ignoring the correc¬ 

tion can induce an uncontrolled systematic error. 

The third method is to combine the first two methods, 
by fitting both the ratios and their sum together (denoted 
as “mixed”), 

-b , (34) 

-Am 

SR^^^{t2,tl) = t2Co -b ^ ^ g_Am ^2) 

+ {t 2 - l)C3e-^’”‘^ + C 4 , (35) 

where Co,1,2,3 and Am are the same as that in the “2- 
state” fit, and C 4 is for the constant contribution from the 
transition between higher excited states and the ground 
state. 


The “2-state” ht makes fully use of the ratios, while it 
is unstable when the position of the plateau is hard to de¬ 
termine (such as for g^). The “sum” fit provides a stable 
estimate of the ground state matrix element, but it suf¬ 
fers from the systematic error from ignoring the 
correction. By combining them together, we can obtain 
a stable fit of all the quantities discussed in this work 
independently, and don’t have to use a joint fit with sev¬ 
eral quantities. The x^/d-o.f. of different quantities and 
quark masses vary between I.O and 1.5 with 18 degrees of 
freedom, corresponding to p-values in the range of [0.08- 
0.46]. The value of Am we obtained at the unitary point 
has a strong dependence on the quantity and varies from 
400 MeV to IGeV. 

The values for the renormalized isovector axial vector 
coupling g\, scalar coupling g^ and quark momentum 
fraction {x)u-d from the three methods at the unitary 
point are listed in Table II. 


TABLE II: Isovector axial-vector coupling g\, scalar coupling 
gg and quark momentum fraction {x }^_4 at the unitary point 
from three fitting methods. See the following three subsec¬ 
tions for the details. 



2 -state 

sum 

mixed 

9a 

1.189(20) 

1.157(18) 

1.166(19) 

9s 

0.61(6) 

0.78(6) 

0.74(4) 

{x}S-d 

0.209(12) 

0.190(13) 

0.193(19) 


A. Vector and Axial vector case 
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FIG. 6: The vector renormalization factor from the charge vs 
the pion mass, from three kinds of fitting methods: 2-state fit 
(red squares), summed slope (black triangles), and the mixed 
fit which combines those two methods (blue dots). The results 
from these different methods are consistent while that from 
the mixed method provides the best signal. 
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The lattice renormalization of the vector current can 
be defined from normalizing the vector charge, 

b ^ Tt[T^{P\ J ^ J_ 

- Tr[r«(P|P)] Zy ^ ^ 

where superscript b is for bare value, and T® = (1+74)72 
is the unpolarized projection operator. Fig. shows 
that all the fitting methods mentioned in the last sec¬ 
tion provide consistent results, while the results from 
the “mixed” method have the best signals among the 
three methods. A constant fit for the cases with ttItt S 
(0.25, 0.4) GeV gives the value of the vector renormaliza¬ 
tion factor as 1.096(6) which is just slightly smaller than 
the value 1.105(4) obtained from the axial Ward identity 

m¬ 



1-12/2 


FIG. 7: The sink-source separation dependence of the matrix 
element of the isovector axial-vector current. The data of the 
isovector case with t 2 = 8 (black squares) are slightly smaller 
than that from the other two separations, while the result 
with t 2 = 10 (blue dots) is consistent with that with t 2 =12 
(red triangles). 


Then the renormalization of the vector current can be 
used to renormalize the axial-vector matrix element with 
polarized projection. 


9a 


9a 


Ei=i,2,3 Tr[r)"(P| / (fx'ilj{x)-/5^ii^{x)\P)] 


3 Tr[r^(P|P)] 


= 9aZv 


Sz^i,2,3 Tr[P(P| / d^xi;{x)'y5'yii^(x)\P)] 
3 Tr[r®(P| J d^x^l){x)"f4tlJ{x)\P)] 


(37) 


the 


where the superscript b/R stands for 
bare/renormalized value respectively and F™ = 
(1 + 74)7275/2 is the polarized projection operator. 

Using gy^ (instead of that from the axial Ward iden¬ 
tity for pion) to renormalize gA as in Eq. (37) could im¬ 
prove the signal of the renormalized gA by ^20% since 
these two matrix elements are correlated. As observed 


in Fig. the sink-source separation dependence for the 
isovector case is mild, while a curve is observable at the 
right side of the plateau due to a larger excited state con¬ 
tribution from the point interpolation field at the sink. 
This is in contrast to the flatter behavior to the left of 
the plateau where the excited-state contribution is ame¬ 
liorated by the smeared source. In Fig. we plot the 
results of the isovector axial-vector coupling g^ from the 
three fitting methods we mentioned. We note that those 
from the “mixed” method are always between those from 
the other two methods, for all the data points in the 
range of ?7i,r G (0.25,0.4) GeV. The values from the three 
methods at the unitary point are listed in Table [Til Sim¬ 
ilar to other lattice calculations at this pion mass (i.e. 
~ 300 MeV), irrespective of which fit is used, the isovec¬ 
tor axial-vector matrix element, g^'^ is ^10% smaller 
than the experimental value 1.2723l23')|45j. 
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FIG. 8: The isovector axial-vector matrix element vs the pion 
mass, from three kinds of fitting method: 2-state fit (red 
squares), summed slope (black triangles), and the mixed fit 
which combines those two methods (the blue dots). The re¬ 
sults from these different methods are consistent while that 
from the mixed method provides the best signal. 


B. Scalar case 

Similarly, the renormalized scalar matrix element with 
the unpolarized projection of the nucleon can be calcu¬ 
lated by. 


ZsTr[F®(P|/ d^xi;{x)^{x)\P)] 
Tr[r'=(P|P)] 


(38) 


where the renormalization constant Zs is obtained from 
the RI/MOM scheme and its value on the ensemble we 
use here is calculated to be 1.1397(54) [41]. On the other 
hand, if one just focuses on the irNa term, 2ZjyimqZsgg, 
the renormalizations of the quark mass Zm and that of 
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both the u and d quark increase as niq decreases, but the 
isovector scalar matrix element is not far from unity over 
the entire quark mass region from light to heavy. This 
has been interpreted to be related to the Gottfried sum 
rule violation |46j where it is found experimentally that 
there are more d antipartons than u antipartons.. 
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FIG. 10: Isovector scalar matrix element vs. pion mass 
from three kinds of fitting method: 2-state fit (red squares), 
summed slope (black triangles), and the mixed fit which com¬ 
bines those two methods (blue dots). The results from these 
different methods are slightly different. 


2-state ^ 

■Q- 

— 

sum ^ 

-A- 


mixed — 

-e- 






[ 


[ 

1 


■ 


■ 


.1 0.2 0.3 0.4 0 

(GeV) 


FIG. 9: The separation dependence of the matrix element of 
scalar current, for both the isovector and the GI part of the 
singlet case. The dependence is mild for the isovector case 
(the upper panel), while obvious for the Cl part of the singlet 
case (the lower panel). 


the scalar matrix element Zs are canceled and so the 
irNa term is free of the renormalization. 

It is interesting to point out that the Cl part of the 
scalar singlet matrix element has a strong sink-source 
separation dependence, as seen in the lower panel of 
Fig.0 At the same time, such a separation dependence 
seems to be canceled between the u and d quarks, so that 
the isovector case in the upper panel of Fig. |^has only a 
mild separation dependence. The results for the isovec¬ 
tor scalar matrix element from the three fitting methods 
are plotted in Fig. [T0|and those at the unitary point are 
listed in Table [HI This shows that, despite the fact that 
there are 2 u valence quarks and only one d quark in the 
proton, the d contribution to the scalar matrix element 
per quark is more than that of the u, as 


9 lci 

‘^9s,ci 


0.67(2) 


(39) 


is much smaller than one. The scalar matrix elements of 


C. Quark momentum fraction 


The quark momentum fraction in the nucleon can be 
calculated with the traceless part of the energy momen¬ 
tum tensor, and it should be consistent between calcula¬ 
tions with two different operators. The first one uses the 
combination of the diagonal temporal and spatial com¬ 
ponents of the energy momentum tensor, 


E ^ Trr(P|/dW(x)|P)] 
^ ^ £;Tr[r®(P|P)] 


(40) 


where 0^{x) = '(/'(a;) 5(74-D 4 - i Ei=i,2.3 7*-D i)i/'(a:) is 

the traceless part of the energy momentum tensor 
and is a measure of the quark fraction of the nucleon 
mass or energy. The related matrix element can be cal¬ 
culated in the rest frame and, as a result, it will have 
a good signal. On the other hand, the operator it¬ 
self can have mixing with lower dimension operators like 
the dimension-3 scalar operator tjj{x)ijj{x). Nevertheless, 
such a mixing will be canceled due the subtraction of the 
diagonal spatial components in Oe- 

The other approach uses the forward off-diagonal ma¬ 
trix components of the energy momentum tensor (r4i) in 
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FIG. 11: The plateau fit values (t 2 = 10 case) of the isovector 
momentum fraction for (a;)E in the rest frame (red square) 
and also (a:)p in a moving frame with different momenta (blue 
dots). The results from both the diagonal and off-diagonal 
components (and also that from different momenta based on 
the off-diagonal matrix components) are consistent, but the 
first approach provides much better SNR. 


a moving frame, 

p ^ Tr[r^(P| /d3a;V:(a;)|(7i^4 + 74^0V^(a^)|J’)] 
“ p,Tr[r-(P|P)] 

(41) 

with Pi being the f-th component of the nucleon momen¬ 
tum. Therefore, it is a measure of the quark momentum 
fraction in a moving nucleon. Such a scheme is free of 
mixing of the lower dimension operators due to its tensor 
structure, while the corresponding matrix element is pro¬ 
portional to the momentum and is thus more noisy than 
that from the first approach, because mixed momentum 
sources are involved for the matrix element of the nucleon 
at non-zero momentum. 

Fig. [TT] shows the plateau fit values of the ^2=10 
case for the quark isovector momentum fraction. They 
are (a;)^ from the diagonal components of the energy- 
momentum tensor with the nucleon in the rest frame and 
also (x)^ from the off-diagonal components in a moving 
frame with different momenta. The results from both 
the diagonal and off-diagonal components (and also those 
from different momenta) are consistent, but (a;)® pro¬ 
vides much better SNR. The sink-source separation de¬ 
pendence is shown in Fig. |12[ for both results based on 
the diagonal components and off-diagonal components. 
It is interesting to observe that the separation depen¬ 
dence of the isovector quark momentum fraction based 
on the off-diagonal components seems to be milder than 
that based on the diagonal ones, for the cases with t 2=8 
and 10. The (a;)^ case with ^2=12 seems to have some 
t dependence at the smeared source end, but it could be 
due to the statistical fluctuation due to relatively poor 
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FIG. 12: The sink-source separation dependence of the isovec¬ 
tor quark momentum fraction for the case of the diago¬ 
nal components of the energy momentum tensor (the upper 
panel) and that of the off-diagonal ones (the lower panel). 


signal. 

As in Ref. m, the renormalization factor for the en¬ 
semble we used has been obtained with the one-loop lat¬ 
tice perturbative theory, as 1.049(3), in the MS scheme 
at 2 GeV. The error is from the uncertainty of the lattice 
spacing. The renormalized values of the isovector quark 
momentum fraction of (x)e from the three fitting meth¬ 
ods are plotted in Fig.[^ and those at the unitary point 
are listed in Table im 


VI. SUMMARY 

We have introduced a new method to calculate the nu¬ 
cleon matrix elements in the connected insertion. The 
stochastic sandwich method (SSM) with low-mode sub¬ 
stitution (LMS) is an approach which uses low modes for 
the all-to-all quark propagator between the current and 
the sink and the corresponding high-mode contribution 
is taken care of by the noise propagator from the sink 





















14 


0.4 
0.35 
0.3 

^ 0.25 

0.2 

^ 0.15 

0.1 
0.05 
0 

0.1 0.2 0.3 0.4 0.5 

(GeV) 


2-State - 

— & 

sum - 

—A- 

mixed - 

—©. 



This shows that, despite the fact that there are 2 u va¬ 
lence quarks and only one d quark in the proton, the d 
contribution to the scalar matrix element per quark is 
more than that of the u, as 


0.67(2) (44) 

‘^9s,ci 

is much smaller than one. This has been interpreted [15] 
to be related to the Gottfried sum rule violation jlTj 
where it is found experimentally that there are more d 
antipartons than u antipartons. 

In the isovector quark momentum fraction case, the 
bare value we obtained at the unitary point on the en¬ 
semble mentioned above is 


FIG. 13: The isovector quark momentum fraction ((x)® in 
the rest frame) vs. the pion mass, from three htting methods: 
2 -state fit (red squares), summed slope (black triangles), and 
the mixed fit which combines those two methods (blue dots). 
The results from the three methods are consistent. 


to the current. We have shown that it is more efficient 
than the sink- and current- sequential methods. How¬ 
ever, it does not scale well with volume which requires 
more low eigenmodes. It will lose its advantage when 
the overhead from calculating the LMS for all the quark 
propagators involved is more than the amount it saves 
compared with the sink-sequential or current-sequential 
method. But this will occur only at volumes much larger 
than that used here. 

We have used three fitting methods. One is a two-state 
fitting including the contamination from the excited-state 
transition and the second is the summed-slope method. 
The third is a mix of these two methods. 

The proton isovector axial-vector coupling we ob¬ 
tain with the overlap fermion at the unitary point with 
771^=330 MeV is 

g\ = 1.166(19) (42) 

which is is ^ 8% smaller than the experimental value. 

The separation dependence of this quantity is mild. 
Since it is smaller than the experimental value on this 
lattice, it is essential to repeat the calculation of g\ on 
larger volumes and with lighter quark masses. 

For the isovector scalar matrix element in the pro¬ 
ton, the renormalized value at MS{2GeV) at the unitary 
point is 

9s = 0.74(4). (43) 


(x),_d = 0.192(19), (45) 

with the renormalization factor 1.049(3) from one-loop 
lattice perturbative theory m- This value is similar to 
those from most lattice calculations laiZKIlHIll and is 
larger than the experimental value. However, the O(a^) 
error has not been considered. It can be assessed by 
imposing the momentum and angular momentum sum 
rules at finite lattice spacing as is demonstrated in a 
quenched calculation [52|. We will return to this issue 
when the complete lattice simulation of the momentum 
and angular-momentum decompositions is carried out. 

We will perform calculations with physical sea quark 
masses in the future. 
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